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ABSTRACT 

We predict the redshift distribution of Gamma-Ray Bursts (GRBs) assuming that they trace the 
cosmic star formation history. We find that a fraction ^ 50% of all GRBs on the sky originate at a 
redshift ^ ^ 5, even though the fraction of the total stellar mass formed by z ~ 5 is only ^ 15%. 
These two fractions are significantly different because they involve different cosmological factors when 
integrating the star formation rate over redshift. Hence, deep observations of transient events, such 
as GRB afterglows or supernovae, provide an ideal strategy for probing the high-redshift universe. We 
caution, however, that existing or planned flux-limited instruments are likely to detect somewhat smaller 
fractions of high redshift bursts. For example, we estimate that the fraction of all bursts with redshifts 
z ^ 5 is ~10% in the case of the BATSE instrument, and ~25% in the case of Swift. We also show that 
the intrinsic distribution of GRB durations is bimodal but signiflcantly narrower and shifted towards 
shorter durations than the observed distribution. 

Subject headings: Cosmology: theory — early universe — gamma rays: bursts 



1. INTRODUCTION 

Gamma-Ray bursts (GRBs) are the brightest electro- 
magnetic explosions in the universe (for a recent review, 
see Piran 2000). Popular models for their central engine 
divide into two main classes: (i) the collapse of a massive 
star to a black hole (BH) (MacFadyen, Woosley, & Heger 
2001, and references therein); (ii) the coalescence of a bi- 
nary system involving a neutron star (NS) and a BH or a 
NS as a companion (e.g. Eichler et al. 1989; Janka et al. 

1999) . The observed association of long-duration GRBs 
with star forming regions (Djorgovski et al. 2001c, and 
references therein) , and the possible supernova signatures 
in rapidly-decaying afterglows (Bloom et al. 1999; Kulka- 
rni et al. 2000; Reichart 2001) favors the first class. Both 
classes of models associate GRB progenitors with compact 
objects (BH or NS) that are the end products in the evolu- 
tion of massive stars. Hence, the GRB formation history 
is expected to follow the cosmic star formation history 
(Totani 1997, 1999; Wijers et al. 1998; Blain & Natara- 
jan 2000) up to the highest redshifts (z ~ 20) at which 
the first generation of stars may have formed (Barkana & 
Loeb 2001). GRBs might therefore provide an ideal probe 
of cosmic star formation at all redshifts that in particular 
is unaffected by dust obscuration (e.g., Blain & Natara- 
jan 2000; Porciani & Madau 2001). In fact, the top-heavy 
initial mass function (IMF) predicted for the first stars 
(Bromm, Coppi, & Larson 1999, 2002; Abel, Bryan, & 
Norman 2000, 2002, Nakamura & Umemura 2001) favors 
massive stars which are the likely source of GRB progeni- 
tors. 

GRB afterglows provide a unique probe of the high red- 
shift universe (Lamb & Reichart 2000; Ciardi & Loeb 

2000) . The bright, early optical- UV luminosity of a GRB 
afterglow is expected to outshine its host galaxy, even more 
so at high redshifts when the typical galaxies are less mas- 
sive than their present-day counterparts (Barkana & Loeb 



2001). The broad-band afterglow spectrum extends into 
the far UV and so the absorption features imprinted on 
it by the intervening intergalactic medium (IGM) can be 
used to infer the evolution of the neutral hydrogen fraction 
and the metal abundance of the IGM during the epoch 
of reionization. In difference from galaxies and quasars, 
which fade rapidly with increasing redshift due to the in- 
crease in their luminosity distance, GRB afterglows main- 
tain an almost constant infrared flux with increasing red- 
shift at a fixed time lag after the GRB trigger in the ob- 
server frame (Ciardi & Loeb 2000). This follows from 
the cosmological time-stretching of the afterglow transient 
(which is intrinsically brighter at earlier times) and from 
a favorable X-correction in the afterglow spectrum. 

The Swift satellite^, planned for launch in 2003, is ex- 
pected to localize roughly one GRB per day. Sorting out 
the subset of all GRBs which originate at high redshifts 
(z ^ 5) would be of particular interest. Observers may em- 
ploy a simple strategy for this purpose. Photometric data 
from a small telescope should be used at first to identify 
those GRBs which possess a Lya trough at a wavelength 
of 0.73/xm(l + z)/6 due to absorption by the IGM. Follow- 
up spectroscopy of those GRBs could then be done on a 
10-m class telescope. In designing this observing strategy 
it is important to forecast which fraction of all GRBs orig- 
inate from different redshifts. For example, it would be 
impractical to search for those very high redshifts which 
amount to a fraction smaller than 10"'^ of all GRBs, be- 
cause barely a single one of them would be found by Swift 
over several years of operation. 

In this paper, we use existing observational and theoret- 
ical work on the cosmic star formation history to predict 
the fraction of all GRBs that are expected to originate at 
different redshifts. In order to keep our results general, we 
make predictions about all GRBs without reference to the 
detection threshold or redshift horizon of any particular 
instrument. To ascertain, however, what the BATSE and 
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Swift instruments are expected to detect, we in addition 
estimate the redshift distributions for these flux-limited 
surveys. 

In §2, we calculate the collapsed fraction of baryons as a 

function of redshift based on the Prcss-Schcchtcr formal- 
ism, and infer the corresponding redshift distribution of 
GRBs. In §3 we use the inferred redshift distribution of 
GRBs to convert the observed distribution of GRB dura- 
tions into the corresponding intrinsic distribution, under 
the simple assumption that its normalized form is redshift 
independent. Finally, we discuss the implications of our 
results in §4. 

2. STRUCTURE FORMATION MODEL 

2.1. Star Formation History 

We adopt the popular view that the formation of cos- 
mic structure has progressed hierarchically from small to 
large scales, according to a variant of the cold dark matter 
(CDM) model. Specifically, we assume a ACDM model 
with density parameters in matter Qm = 1 — il\ = 
0.3 and in baryons fie = 0.045; a Hubble constant of 
h = i?o/(100kms~^Mpc~^) = 0.7, and a scale-invariant 
power spectrum of density fluctuations with an amplitude 
(Tg = 0.9 on a scale of 8/i~^Mpc. 
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Fig. 1. — History of cosmic comoving star formation rate (SFR) in 
units of M0yr~^Mpc~^ as a function of redshift. Solid line: cool- 
ing due to atomic hydrogen only; Dashed line: added cooling via 
molecular hydrogen. The star formation efficiency is assumed to be 
V* = 10%. 

Our star formation model closely follows that of Barkana 
& Loeb (2000), and here we only briefly describe the key 
assumptions (see also Santos, Bromm, & Kamionkowski 
2002). The abundance and merger history of the CDM 
halos is described by the extended Press-Schechter formal- 
ism (Lacey & Cole 1993). We assume that the IGM has a 
two-phase structure, consisting of a neutral and an ionized 
phase. The fraction of the cosmic volume flUed with H II 
regions is given by 
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exp[-0.495(.2 - 5.6)] 



for z > 5.6 
otherwise 



(1) 



This assumed ionization history fits the semi-analytical 
calculation of Barkana & Loeb (2001) and is consistent 
with numerical simulations of reionization (Gnedin 2000, 
2001; Razoumov et al. 2001) and the latest data on 
quasars in the redshift interval 5 <J z < 6.3 (Becker et 
al. 2001; Djorgovski et al. 2001a; Fan et al. 2001). At 
high redshifts, z ^ 20, the universe is predominantly neu- 
tral. Once the first luminous objects form, an increasing 
fraction of the IGM becomes ionized. At Zrcion ~ 7, the 
ionized phase in our model comprises a volume fraction 
of ~ 50%, and reionization of the IGM is complete by 
z w 5.6. 

Within each phase of the IGM, stars arc able to form 
in two different ways. The first mechanism pertains to 
primordial, metal-free, gas. Such gas undergoes star for- 
mation provided that it falls into a sufficiently deep CDM 
potential well, or equivalently, into a CDM halo more mas- 
sive than a critical value. For the neutral medium, this 
minimum mass is set by the requirement for the gas to 
cool. Radiative cooling by molecular hydrogen (H2) al- 
lows star formation in halos with a virial temperature 
Tvii ^ 300 K, while atomic cooling dominates for halos 
with Tvir ^ 10^'^ K. The corresponding minimum circular 
velocities are Vc ^ 2.5 km s~^ and ~ 12 km s~^, respec- 
tively. Since H2 can be easily photo-dissociated by pho- 
tons below the Lyman- limit, its significance in the cosmic 
star formation history is unclear (e.g. Haiman, Abel, & 
Rees 2000; Ricotti, Gnedin, & Shull 2001, and references 
therein), and so we show results with and without H2 cool- 
ing. These two theoretical models are likely to provide 
conservative bounds for the true star formation history at 
^; 5. The construction of more tightly constrained mod- 
els has to await further advances in our understanding, 
both observational and theoretical, of star formation at 
the highest redshifts. 

For the ionized medium, on the other hand, the mini- 
mum threshold mass is given by the Jeans mass, since the 
infall of gas and the subsequent formation of stars requires 
that the gravitational force of the collapsing CDM halo be 
greater than the opposing pressure force on the gas. Af- 
ter reionization, the IGM is photo-heated to temperatures 
^ lO'^K, leading to a dramatic increase in the Jeans mass. 
We model the suppression of gas infall according to results 
from spherically-symmetric collapse simulations (Thoul & 
Weinberg 1996). Expressing the Jeans mass as the equiva- 
lent halo circular velocity, we assume complete suppression 
for halos with Vc < 35 km s"-'^, no suppression for Vc > 93 
km s~^, and a linear interpolation in between so that ~ 
50% suppression occurs at Vc ^ 55 km s~^. 

Within our model, the second mechanism to form stars 
occurs in gas that has experienced a previous burst of 
star formation, and is therefore already somewhat enriched 
with heavy elements. Such gas, residing in a halo of mass 
Ml, can undergo induced star formation triggered by a 
merger with a sufficiently massive companion halo of mass 
M2 > 0.5Mi. We finally assume that stars form with an 
efficiency of ry, ~ 10%, independent of redshift and regard- 
less of whether the gas is primordial or pre-enriched. This 
efficiency yields roughly the correct fraction of Qb found 
in stars in the present-day universe. 

Figure 1 shows the resulting star formation histories. 
Our theoretical models agree well with observational esti- 
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mates of the cosmic star formation rate (SFR) at z <J 2 
(e.g., Blain et al. 1999). It is evident that there are two 
distinct epochs of cosmic star formation, one at 3, 
and a second one at 2; ~ 8 (or at even higher redshifts 
if H2 coohng is effective). Again, we emphasize that the 
true history of the cosmic SFR is Ukely to he between the 
two curves in Figure 1, depending on how complete the 
destruction of H2 as a function of redshift is. 

In deriving the redshift distribution of GRBs in the next 
section, we do not make any assumptions on the possible 
variation of the IMF for stars forming at different red- 
shifts. Instead, we only assume that baryons are incorpo- 
rated into stars, regardless of their specific properties, with 
the overall rate calculated in this section. Let us, however, 
briefly discuss the possible differences in star formation at 
high and low redshifts, based on recent theoretical work 
implying that star formation at high redshifts might have 
proceeded very differently from the present-day case, lead- 
ing to stars with typical masses of M* ^ lOOM© (Bromm 
et al. 1999, 2002; Abel et al. 2000, 2002; Nakamura & 
Umemura 2001). After the first stars have formed, the 
subsequent generation of stars forms out of gas that was 
already enriched with heavy elements. This enriched gas 
could have cooled more efficiently, and was able to reach 
lower temperatures. Star formation, then, is expected to 
result in a less top-heavy IMF. As shown by Bromm et 
al. (2001), the transition from a top-heavy to the more 
standard (Salpeter) IMF occurs when the mass fraction 
in metals exceeds a critical value of ~ 10~^Zq. Gas with 
a metal abundance below this threshold is therefore still 
expected to form very massive stars. An IGM metal abrm- 
dance of ^ 10^"^ Zq approximately corresponds to the pro- 
duction of enough ionizing stellar photons to reionize the 
universe. Star formation ad, z '^7 might consequently have 
been dominated by very massive stars, whereas at lower 
redshifts, stars form with an IMF close to the Salpeter 
form. 

2.2. Redshift Distribution of GRBs 

Under the assumption that the formation of GRBs fol- 
lows closely the cosmic star formation history with no 
cosmologically-significant time delay, we write for the 
number of GRB events per comoving volume per time: 
■0grb(-z) = Vgrb X i'*iz), where ip*{z) is the stellar mass 
produced on average per comoving volume per time, as 
calculated in §2.1. The efficiency factor, t/grb, links the 
formation of stars to that of GRBs, and is in principle a 
function of redshift as well as the properties of the underly- 
ing stellar population. Massive stars of Population I differ 
fundamentally from those of Population III; moreover, it 
is at present not well understood how a massive Popula- 
tion III star may give rise to a GRB (see Fryer, Woosley, 
& Heger 2001; Schneider et al. 2001). Given the current 
state of uncertainty with regard to the central engine of 
GRBs (e.g., Piran 2000), we make the simplifying assump- 
tion that both populations of massive stars are connected 
to GRBs in a similar way, and take tjgrb to be indepen- 
dent of redshift. While this simplifying assumption follows 
from the lack of better information, our analysis provides 
a starting point for future improvements as soon as better 
observational constraints on high redshift GRB and star 
formation will become available. 

If we now consider a time interval Atobs in the observer 



frame, the total number of GRBs, regardless of whether 
they are actually observed or missed, can be written as 



N{> z) 



^GRB(.)^^^d. , (2) 



where dV/ dz is the comoving volume element per unit red- 
shift, given by 
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The luminosity distance, d^, to a source at redshift z is 

dt 

(1 + 2') 

with 



dL = C(l + Z) [ (1 + z') 
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dz 



{i + z)HoVnm{i + zf + nA (5) 



in a flat universe. The fraction of bursts that originate 
from a redshift of z or higher. /(> z) = N{> z)/N{> 0), 
is independent of constant parameters such as Aiobs, Vgrb 
or the beaming factor of the GRB emission. The integrand 
in equation (2) contains the differential comoving volume 
element, dV/dz, as is appropriate for the calculation of 
an event rate. Since GRB events are communicated via 
photons, we integrate over the redshift-dependent comov- 
ing volume element along our past light cone. We observe 
these events over a fixed time window, Aiobs, which cor- 
responds to Aiobs/(l + z) in the source frame. If, on the 
other hand, we were interested in determining the amount 
of stellar fossils that have accumulated over cosmic time 
in a local comoving volume (see below), we would have to 
simply integrate over cosmic time along our past worldline. 
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Fig. 2. — Redshift distribution of GRBs. (a) Fraction of bursts 
that originate at a redshift higher than 2 vs. 2. The curves cor- 
respond to the two star formation histories in Figure 1. The data 
points reflect ^ 20 observed redshifts (from GhiseUini 2001). (5) 
Fraction of bursts per logaritmic interval of (1 + 2) vs. 2. The 
curves have the same meaning as in panel (a). Dotted lines: Red- 
shift distribution of the baryonic mass fraction in stars, calculated 
for the case of atomic hydrogen cooling. 
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Fig. 3. — Redshift distribution of all GRBs as compared to that 
measured by flux-limited surveys. The convention for the panels 
and the meaning of the symbols is the same as in Fig. 2. Solid 
lines: All GRBs for star formation through atomic line cooling. Dot- 
dashed lines: Expected distribution for Swift. Long- dashed lines: 
Expected distribution for BATSE. Note that the curves for the two 
flux-limited surveys are very uncertain due to the poorly-determined 
ORB LF. 

In Figure 2, we show /(> z) together with the differen- 
tial distribution (d//dz) for the two star formation histo- 
ries of Figure 1 . It is evident that a significant fraction of 
all bursts is predicted to occur at high redshifts, namely 
f{z > 5) ^ 50%, and that a few percent of all bursts occur 
at redshifts as high as z ^ 20. Evaluating the mean red- 
shift for GRBs using the distributions of Figure 2, we find 
^ ~ 5 in the case of atomic cooling and z ^ 8 for molecular 
cooling. Currently, only ^ 20 bursts have known redshifts 
(Kulkarni et al. 2000; Djorgovski et al. 2001c; Ghisellini 
2001), and we include this small observed sample in Figure 
2. The fact that the data points lie below our theoretical 
prediction could be due to small-number statistics as well 
as due to a redshift-dependent incompleteness bias. 

We stress that equation (2) gives the fraction of tran- 
sient events observed on the sky, and not the fraction of 
all baryons that have been incorporated into stars by a 
redshift z relative to the same fraction today. This latter 
quantity is given by 



/*(> z) 



J^Mz')\dt/dz'\dz' 
J^Mz')\dt/dz'\dz' ' 



(6) 



and is shown in Figure 2 for the case of atomic cooling. 
As can be seen, the fraction of all stars that are formed 
at 2; ^ 5 is ~ 15%. Collecting photons from our past light 
cone is therefore a more efficient way of probing the first 
stars than sorting through the fossil stellar record in the 
present-day universe. 

It is important to emphasize that the analysis presented 
here pertains to all bursts, regardless of whether existing 
or previous instruments have actually been able to detect 
them. If the horizon of previous instruments was limited 
to 2 <C 5, then our predictions provide important moti- 
vation for the construction of more-sensitive instruments 
that would trigger on GRBs out to the highest redshifts. 



The fraction of all bursts that were detected by any given 
instrument depends on the instrument-specific flux sensi- 
tivity threshold and on the poorly-determined luminosity 
function (LF) of GRBs (see e.g., Schaefer, Deng, & Band 
2001; Schmidt 2001; Norris 2002). 

It is nevertheless instructive to ascertain what existing 
or planned instruments like BATSE and Swift are expected 
to find. To this extent, we modify the GRB event rate to 
have 



/>oo 

V'GRb(2;) = ?7GRB'0*(^) / p{L)dL 



(7) 



Here, p{L) is the GRB LF with L being the intrinsic pho- 
ton luminosity (in units of photons s~^). If fum denotes 
the sensitivity threshold of a given instrument (in photons 
s" 



cm ), then the minimum luminosity is 



(8) 



This expression is derived with a spectral index of a = 2 
for L oc (Band et al. 1993). For definiteness, we as- 
sume a log- normal distribution function (e.g., Woods & 
Loeb 1995) 



^^exp[-(ln(L/Lo))V(2a^)]^ > (9) 



where a and Lq are the width and the average luminos- 
ity, respectively. Recently, Sethi & Bhargavi (2001) have 
shown that both the observed number count-flux rela- 
tion as well as the existing afterglow redshift data are 
consistent with a log-normal LF for best-fit parameters 
(taking into account the effect of beaming): a = 2 and 
Lq = 2 X lO^^s"^, and we adopt these values in the follow- 
ing analysis. To determine the expected redshift distribu- 
tion as observed by BATSE and Swift, we use equation (2) 
together with the GRB rate given in equation (7). The fiux 
thresholds are /lim = 0.2 and 0.04 photons s~^ cm^^ for 
BATSE and Swift, respectively (Lamb & Reichart 2000, 
and references therein). In Figure 3, we show the same 
quantities as in Figure 2, but now comparing the distribu- 
tions for BATSE and Swift with our theoretical prediction 
for atomic line cooling. It can be seen that in the case 
of BATSE a fraction of f{z > 5) ^ 10% of all bursts 
originates from high redshifts, whereas the corresponding 
fraction for Swift is f{z > 5) ^ 25%. We emphasize again 
that these numbers are uncertain due to the poorly-known 
GRB LF. Figure 3 nicely demonstrates the asymptotic 
character of our theoretical prediction, pertaining to a fu- 
ture 'ultimately-sensitive' instrument. Indeed, using the 
LF above, we estimate that an instrument with a sensitiv- 
ity of ~ 50 times better than Swift would be able to detect 
the full theoretically-possible sample of bursts from z ^5. 

The detcctability of z ^ 5 GRBs is also a crucial ingre- 
dient in estimating the fraction of all well-localized bursts 
that have no detectable optical afterglow, the so-called 
"dark GRBs". Various authors have used the fraction 
of dark bursts in the currently observed sample of GRBs 
to constrain the amount of dust obscured star formation 
(e.g., Djorgovski et al. 2001b). The resulting fraction of 
dark GRBs estimated for different redshifts depends on 
the, presently unknown, level of incompleteness in the ob- 
served sample. In the context of our model, we predict 



5 



that all GRB afterglows originating at z > 6 arc optically 
dark. The intervening, partially neutral IGM would cfB- 
ciently absorb the rest-frame UV afterglow that would oth- 
erwise have been redshifted into the optical band (see also 
Fruchter 1999; Piro et al. 2002). These bursts might give 
rise to the recently discovered class of X-ray rich GRBs 
(e.g., Piro et al. 2002; see also Schneider et al. 2001) due 
to the redshifting of the source- frame 7-rays into the X-ray 
band. 

3. BURST DURATIONS 

The duration of a GRB reflects the characteristic 

timescale over which the central engine is active and is 
therefore a diagnostic of the GRB progenitor. The dis- 
tribution of GRB durations has been determined by the 
BATSE instrument on board the Compton Gamma Ray 
Observatory (as summarized in Paciesas et al. 1999), 
and is observed to be bimodal with a population of short 
bursts centered on Tobs ^ 0.3 s, and long bursts around 
^obs ~ 30s (Kouveliotou et al. 1993). For the definition 
of the burst duration, Tobs, we use the interval of time 
over which a GRB contains from 5 to 95% of its total ob- 
served 7-ray counts, also denoted as Tgo in the literature. 
Since bursts originate over a broad range of redshifts, the 
question arises as to what the intrinsic distribution of du- 
rations is like. For simplicity, we assume in this section 
that BATSE was capable of sampling the full redshift dis- 
tribution of GRBs shown in Figure 2. This provides us 
with the maximum level of distortion that cosmological 
time dilation could have had on the observed distribution 
of burst durations. The true intrinsic distribution of dura- 
tions for the BATSE-triggered bursts should lie in between 
the observed distribution and the intrinsic one calculated 
in this section. 

The number of observed bursts in a given bin i, with 
an observed duration Tobs i 

and a width ATobsi, can be 

written as 

/■°° dP 1 df 

AATobs.i = A^tot ATobs,i / :T?f^ (T) -— - -f iz)dz , 
Jo al (1 -I- z) az 

(10) 

where A^tot = X^i AAobs.i is the total number of bursts 
in the sample. We assume that the intrinsic distri- 
bution, dP/dT, is independent of redshift and satisfies 
/p°°(dP/dT)dT = 1 and (dP/dT) > 0. The intrinsic burst 
duration, T, is related to the observed one by the cosmo- 
logical time dilation, T = Tobs,i/(l + z), and (df/dz) is 
the GRB redshift distribution, as calculated in §2. We 
replace the integration in equation (10) by a summation, 
covering the range of intrinsic durations, T, with the same 
number of bins, Abin = 23, as the observcxl histogram. 
The inversion problem is then uniquely defined. We carry 
out the deconvolution with the standard iterative Lucy 
method. This is a reliable technique, derived from Bayes' 
theorem, to solve a set of linear equations with additional 
constraints on the unknowns (Lucy 1974). The stability 
of the algorithm is improved by limiting the change in the 
unknowns in each iteration, and by smoothing over adja- 
cent bins. To this extent, we use equation (lib) in Baugh 
& Efstathiou (1993) with parameter values oi (3 = 0.8 and 
e = 0.9. We have verified that the solutions are not very 
sensitive to the choice of these parameters. 

The result of this inversion is shown in Figure 4, where 



we compare the derived intrinsic distribution to the ob- 
served one, dP/drobs(T'obs,i) = AAobs,j/(A"tot APobs.i)- It 
is evident that the intrinsic durations are systematically 
shifted to shorter values due to the cosmological time di- 
lation. The bimodality is preserved, with peaks that are 
narrower than the observed ones (note that the horizon- 
tal scale is logarithmic). The two star formation histories 
discussed in §2 lead to similar intrinsic distributions. The 
shift to shorter durations, however, is more pronounced in 
the case of star formation via H2 cooling. The mean in- 
trinsic durations characterising the first, short-duration, 
peak are ~ 0.05 s for cooling due to atomic hydrogen, 
and ~ 0.03 s for H2 cooling. The corresponding num- 
bers for the long-duration peak are ~ 7s and ~ 5s, re- 
spectively. These differences in the mean durations arc a 
direct consequence of the fact that GRBs originate on av- 
erage at somewhat higher redshift if H2 cooling is effective. 
Note that, statistically, the longest duration bursts, with 
Pobs > 1000 s, are expected to originate at high z, and this 
could bo a successful selection strategy for observations 
targeting high-redshift GRBs. The shift to longer dura- 
tions due to the cosmological time dilation could in part be 
compensated by the following subtle selection effect which 
we ignore in this paper. Sources at high z will on average 
have lower fluxes, and observations with a given sensitivity 
threshold will therefore only detect the brightest portion 
of the total emission, thus systematically underestimating 
the true duration of the burst. 
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Fig. 4. — The maximal effect that cosmological time-dilation may 
have on the observed distribution of GRB durations. Distributions 
are shown as probability per logarithmic duration interval vs. du- 
ration (in s). Solid line: Case of atomic hydrogen cooling. Dashed 
line: Case with molecular cooling included as well. The data points 
correspond to the observed distribution by BATSE (Paciesas et al. 
1999). The intrinsic durations are systematically shifted to smaller 
values due to cosmological time dilation. 

Based on observations of GRB afterglows with known 
redshifts. Frail et al. (2001) have recently presented evi- 
dence for a standard amount of energy release in GRBs, 
E ~ 10^^ erg. Making the simplest assumption of a con- 
stant energy for all long-duration bursts (which are the 
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only ones with measured redshifts so far), one can easily 
derive the luminosity function from the intrinsic distribu- 
tion of burst durations. The luminosity of a burst is then 
simply « E/T^ and the resulting luminosity function is 
obtained by inverting the horizontal axis in Figure 4 and 
changing T to E /T . The long-duration bursts would then 
narrowly cluster around a luminosity of ^ lO^^erg s~^. 

4. DISCUSSION 

We have derived the redshift distribution of GRBs out 
to ^ 20 under the assumption that the GRB event rate 
traces the cosmic star formation rate. We find that > 50% 
of all GRBs on the sky originate from a redshift of 5 or 
higher. On the other hand, the fraction of baryons that 
have been incorporated into stars by z '--^ 5 is much smaller, 
comprising only ^ 15% of the stellar mass formed by to- 
day. The difference between the two fractions follows from 
the different cosmological factors in the redshift integra- 
tions for the statistics of transient events on the sky as 
compared to the census of fossil objects in the local uni- 
verse. The favorable statistical bias towards high-redshift 
events on the sky is expected to apply also to Type II su- 
pernova explosions which are related to the formation of 
massive stars in a similar way as GRBs. Despite their dim- 
ming with increasing redshift, high-redshift supernovae 
will be detectable with sufficiently sensitive telescopes such 
as the Next Generation Space Telescope^ (NGST; Miralda- 
Escude & Rees 1997; Woods & Loeb 1998). In fact, our 
calculation implies that without any additional bias (such 
as redshift-dcpendent dust extinction) approximately half 
of all Type II supernovae detected by NGST will originate 
at z ^ 5. Deep observations of high-redshift GRBs and 



supernovae offer an ideal window into the earliest epoch of 
cosmic structure formation. The lengthening of the dura- 
tion of these transients by a factor (1 -f z) makes it easier 
for observers to monitor their lightcurves. 

Different instruments may find GRBs up to different 
redshifts, depending on their detection sensitivity and the 
highly uncertain GRB luminosity function (Schaefer et al. 
2001; Schmidt 2001; Norris 2002). A trigger-unbiased way 
to infer the redshift evolution of the GRB event rate is to 
compare the number counts of GRBs with the same ab- 
solute (intrinsic) luminosity in different redshift bins. If 
future observations of this type were to determine a mean 
redshift for the GRB distribution significantly lower than 
the one predicted in this paper, then this would indicate 
either that GRB formation at high z is substantially sup- 
pressed, or that GRBs originate from the coalescence of 
binaries with a time delay of a few Gyr between the for- 
mation of a massive star and the GRB event. 

Recent observations indicate that a large fraction, ~ 
50%, of all well-localized GRBs have no associated optical 
afterglow, and are classified as "(optically) dark GRBs" 
(e.g., Piro et al. 2002). According to our model, a sub- 
stantial fraction of these dark bursts could originate from 
2^6. The intervening, partially neutral IGM would ef- 
ficiently absorb the rest-frame UV afterglow that would 
otherwise have been redshifted into the optical band. 
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